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We consider diagonal cylindrically symmetric metrics, with an interior representing a general non- 
rotating fluid with anisotropic pressures. An exterior vacuum Einstein-Rosen spacetime is matched 
to this using Darmois matching conditions. We show that the matching conditions can be explicitly 
solved for the boundary values of metric components and their derivatives, either for the interior or 
exterior. Specializing to shearfree interiors, a static exterior can only be matched to a static interior, 
and the evolution in the non-static case is found to be given in general by an elliptic function of 
time. For a collapsing shearfree isotropic fluid, only a Robertson- Walker dust interior is possible, 
and we show that all such cases were included in Cocke's discussion. For these metrics, Nolan and 
Nolan have shown that the matching breaks down before collapse is complete and Tod and Mena 
have shown that the spacetime is not asymptotically flat in the sense of Berger et al. The issues 
about energy that then arise are revisited and it is shown that the exterior is not in an intrinsic 
gravitational or superenergy radiative state at the boundary. 
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I. INTRODUCTION AND SUMMARY 

Many papers have considered cylindrical solutions, 
with or without matching: see e.g. P, [3, 0] and 
Stephani et al. Chapter 22. Here we take diagonal 
metrics in both the interior and exterior, and initially 
allow an anisotropic fluid interior. The exterior is a vac- 
uum Einstein- Rosen (ER) solution 6]. Our aim was to 
develop a solution where one could explicitly see the rela- 
tion between source motion and gravitational radiation, 
albeit in a physically unrealistic case. We were therefore 
interested in collapse, though for non-static cases one can 
easily reverse the sense of time so that collapse becomes 
expansion and vice versa. 

We set out the metrics and the Darmois matching con- 
ditions (which preclude surface shells in the boundary) 
for a timelike boundary in sections |TT] and IIIIl and show 
that the junction conditions can be explictly solved for 
the boundary values for the exterior in terms of interior 
quantities and vice versa. This extends the work of 

Specializing to the shearfree case in section HVl we are 
able to give a first-order ordinary differential equation 
for the time evolution of the interior whose solution is 
in general an elliptic function. It is shown that a static 
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exterior implies a static interior. 

Further specializing to a isotropic fluid in section |Vl 
we prove that only a Robertson- Walker (RW) dust inte- 
rior is possible and that all such interiors are included 
in the discussion of Cocke Q . The matching then leads 
in section IVII to specific behaviour of the ER functions 
at the boundary. However, previous work of Nolan and 
Nolan [9] shows that the matching breaks down before 
collapse is complete and Tod and Mena fTo*! showed that 
the solutions cannot be asymptotically flat in the sense 
of Berger et al. [ll[ , which makes them unsatisfactory for 
our purposes. 

Finally in section IVHI we consider whether there are 
waves in the exterior by asking if there is energy trans- 
port. Because a cylindrically symmetric spacetime can- 
not be asymptotically flat, we cannot employ the usual 
global definition of radiation for isolated bodies due to 
PJl, and various alternatives are discussed. We conclude 
that the exterior of a cylindrical region of a collapsing 
RW dust solution cannot be in an intrinsic gravitational 
or superenergy radiative state at the boundary, and infer 
that no radiation is transferred to or from the interior. 

The specializations made here did not lead to solutions 
of the type we hoped for. This is a consequence of the ad- 
ditional restrictions imposed in the hope of avoiding the 
full complexity of the problem in the general anisotropic 
case. Nevertheless we believe that it is necessary to have 
the results obtained in these more restricted cases as a 
first step: we hope in the future to undertake further 
study of the shearing dust and shearfree anisotropic fluid 
interiors, which may establish whether such results as the 
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breakdown of the matching in the FRW case hold more 
generally. 



II. COLLAPSING ANISOTROPIC FLUID 
CYLINDERS 

We consider a collapsing cylinder filled with 
anisotropic non dissipative fluid bounded by a timelike 
cylindrical surface T, and with energy momentum tensor 
given by 

= (M + Pr)VaVi3+ Prga(3 

+ {P, - Pr)S^Sp + (P^ - Pr)K^Kp, (1) 

where is the energy density, Pj. , and P^ are the prin- 
cipal stresses and Vq, Sa and Ka are vectors satisfying 



V'^Va = -1, S'^Sa = K°'Ka = 1, 



(2) 



We assume the general time dependent diagonal non ro- 
tating cylindrically symmetric metric 



A^idt^ -dr^)+B^dz^ + C^d(l)\^ (3) 



where A, B and C are functions of t and r. To represent 
cylindrical symmetry, we impose the following ranges on 
the coordinates 

-oo<t<oo, < r, -GO < z < cx), < < 27r, (4) 

where we assume C = at r = which is a non-singular 
axis. We number the coordinates x'^ — t, — r, = z 
and x'^ = 4> and we choose the fluid to be comoving in 
this coordinate system; hence from and ^ 



Vo,^~A5l, So^^BSl, Ko^^CSl 



(5) 



Calculating the motion of the fluid according to its ex- 
pansion and shear (Jap, 



(6) 
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= l^(a;/3) + V(^,^V^Vp) - - e{gc,f3 + V^Vp), (7) 



by using Q and (jl]) we obtain for the expansion, 

^ I (a b c\ 

^^A^A + B + cj' 
and for the non zero components of the shear, 
A 
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where the over dot stands for differentiation with respect 
to t. The Einstein fleld equations. Gap = uTap, for ([T]), 
and ([5]) have the non zero components. 



Gn 









BC 


B" 


G" 


A 


1^ 




^ BC 




^ ~C 



^01 — ^ 



A' 


(B' 


G'\ 


B' 


G' 


A 


[b' 


g) 


B 


G 


B' 


C' 


A 


fB' 


G' 


B 


G 


^A 




G 



KflA^ 



B G\A' „ 

b + c)a='^ 



(12) 
(13) 
(14) 



^11 




A' 




G'\ 


B' 


G' 


A 




'g) 


+ B 


G 



BC 

^ BC 

= KPrA^ 



(15) 



G 



22 




A^ C^ 

~A^~C 



= kP,B^ 
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A^ 
~A 



B^ 



= kP^G^, 



(17) 



where the prime stands for differentiation with respect 
to r. 

From ll]), ^ and ([5]) we have the non trivial compo- 
nents of the Bianchi identities T°'p;a = 0, 

+ Pr)j + {p + + (/i + = 0,(18) 



P; + (M + Pr 



A' 



{Pr~P.)^ + {Pr-P^)^=Q.m 



For the exterior we take an ER spacetime Q 

ds\ = ~e^^^-^''>idT^^dR^)+e^-^dz^+e-^-^R^dc^^, (20) 

where 7 and ip are functions of T and R, and for the 
vacuum field equations we have 



and 



V,TT - ■V,R.R -fT = 0, 

H 



(21) 



(22) 



Equation (|2ip is the cylindrically symmetric wave equa- 
tion in an Euclidean spacetime, suggesting the presence 
of a gravitational wave field. 
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III. JUNCTION CONDITIONS 

As the boundary must be comoving with the fluid inte- 
rior, it will be timelike, and given by r ^constant in the 
interior metric ([3]) and a curve R{T) in the ER metric 
(1^0)) . Matching the collapsing cylinder at S to the ER 
spacetime, Darmois' junction conditions 0, give us, 
after a little algebra. 



interior field equations. 

We note also that in order for the surface to be timelike 
we require {dR/dTY < 1. 

One can reorganize these equations, together with (|22p . 
to conversely give the exterior functions' values on the 
boundary in terms of the interior. Rewriting (j23p . and 
differentiating ([M)) and (PS)) with respect to r in the sur- 
face, we obtain 



dR 
dT 



-,1/2 



dT ^ Adt = dT, 



R = BC, 



s {Bcy 

A ■ 



(23) 

(24) 
(25) 
(26) 

(27) 
(28) 



In this form it is easy to see that if the exterior is known, 
the boundary values of B and C, and 1/A times their 
first derivatives with respect to t and r (equivalently, the 
derivatives with respect to proper time in the surface and 
proper distance orthogonal to it), can be solved for. The 
value of A is not fixed, as one can redefine the t and 
r coordinates, but its evolution can be found from the 



B 



e^{Tr^,T + R.r^.B) 



R.r 



~ A' 
s {BC) 
A ' 



(29) 
(30) 
(31) 



Solving pO|) and (p7)) for ip,T and tp,u, and substituting 
in (121), the results are IMl) and 



s B,t{BC),r~B,,{BC),t 
B[{BC)l~{BC)l] ■ 

s B,r{BC),r-B,t{BC),t 

B[{BC)l-{BC)l] ■ 
E AB 
- [(i3C)2 -(BC)2,]V2' 



(32) 
(33) 
(34) 



7,T 



1,R 



S 2 C[B,t{BC),r - B^r{BC)MBC),rB.r - {BC),tB^t] 
B[{BC)^, - (BC)2 ]2 

S C{[B^t{BC),r - BABC),t? + [{BC),rB,r - (BC),tB 



B[{BC)l-{BC)l 



2 12 



(35) 
(36) 



The radius TZ of the collapsing cylinder as measured by 
the circumference in the exterior ER spacetime is given 

by 



(37) 



IV. SHEARFREE COLLAPSING SOLUTION 

Now let the motion of the collapsing cylindrical fiuid 
be shearfree, CTq/j = 0. Then we can integrate (|9] fTT|) and 
obtain 



B = b{r)A, C^c{r)A, 



(38) 



where b{r) and c(r) are arbitrary functions of r and the 
metric ^ becomes 

dsl = A^i-dt^ + dr^ + b^dz^ + c^d(t?). (39) 



We observe that by contrast the shearfree condition in a 
collapsing spherical distribution of matter in a comoving 
frame leaves two unknown functions of time and radial 
coordinate in the metric [l4| . 

Substituting ([M]) into ([T^ we obtain 



A' A A' 

a-'aa='^ 

which can be integrated producing 
A^ ' 



w{t) + a(r) ' 



(40) 



(41) 



where w is an arbitrary function of t and a is an arbi- 
trary function of r. The expansion ([8]) for ([39)1 and ((4T] 

becomes, 



(42) 



which, as for the shearfree isotropic fiuid spherical col- 
lapse, depends also only on t 1^. 



4 



The field equations (|12j|17p with ^ and (gT]) become 



fb' c! 

= 3w + 2(w + a)a" — 3a' + 2(w + a)a — H 

\ h c 



b' c' 



b c be 
nPr = 2(w + a)w - 'dw^ + 3a'^ 
, (b' c' 



(43) 



-2(w + a)a' - + - +(w + a)V-, (44) 
^ b c I be 



kP^ = 2{w + a)w - 3w^ - 2{w + a)a" + 3a'^ 

c' c" 
-2(w + a)a'— + {w + af — , 
c c 

kP^ = 2{w + a)w - 3w'^ - 2{w + a)a" + Sa'^ 

b' b" 
-2{w + a)a'- + {w + a)^y ■ 



(45) 



(46) 



From the junction condition (j26p we have that on the 
boundary (|44p is given by 



2VtVt-itf + 02^"^ +ci^l + co^Q, (47) 



where 



Si = w + a, Co = 3a , ci — —la — H , C2 — -, — ■ 

\ b c I be 



We can integrate (j47|) producing 



(48) 



(49) 



where wq is an integration constant. Thus in general the 
time dependence is given by an elliptic function whose 
parameters are fixed by values at the boundary. 

If the exterior spacetime is static, tp^T = 0, then the 
field equations (|2ip and reduce to the static Levi- 
Civita spacetime, 



and 



(50) 



(51) 



where ai and 02 are integration constants: for a positive 
mass source we need a2 > [l^, [l^. For the solution 
([39| and (|4T|l we have the following relations on S. From 
(|24l) and jSO]) we have 



ai s & 



and from (jM]) and (pS)) 



s be 

R = 



[w + a)2 
With dUl) and ([521) we obtain 

ai(w + a)2°^+i = 6°^+ic°^ 



(52) 



(53) 



(54) 



since w is a function of t this relation can hold only if w 
is constant. Hence, for shearfree cylindrically symmetric 
anisotropic fluids if the exterior spacetime is static, i.e. 
the Levi-Civita spacetime, the cylindrical source must be 
static too. 

Not all the junction conditions have been used, since no 
specific model of the interior has been given (and static 
shearfree fiuid solutions may not all contain a suitable 
matching surface, for instance they may not contain a 
surface where P,. = 0). However, the remaining junction 
conditions can be satisfied for anisotropic and isotropic 
fluids, shells, and other suitable choices of interior, see, 
for example, [H, [13, [H, [H - 



V. CYLINDRICALLY COLLAPSING 
ISOTROPIC FLUID 

If the collapsing cylinder is filled with isotropic fiuid 
then Pr = Pz = P4, and from (|44ll46p we have 



{w + a) 
{w + a) 



„ ,b'\ , Jc" b'cJ 



(9^5) 
= (56) 



Then from (|55p and (|56p . assuming w has non- trivial time 

dependence. 



(57) 



,b' ,c' b" c" b' c' 
a =a — = a—, — = — = — — , 
c b c be 

which can be shown by direct calculation to reduce the 
Weyl tensor to CafSjS — 0, i.e. the spacetime inside the 
cylinder is conformally flat. All conformally flat perfect 
fluid solutions are known, and all are shearfree [5j]. If a 
conformally flat perfect fluid has a barotropic equation 
of state, then (Triimper, cited in [2^) it is RW. We now 
show directly that this is the case, i.e. that the interior 
must be RW and thus conformally fiat, without assuming 
barotropy, using regularity at the axis instead. 
If a' 0, then we must have 

from which it easily follows, using (|57|) . that a, b and 
c are all proportional to e^"^ for some non-zero constant 
[3. This is not consistent with having a non-singular axis 
where C = at r = 0. 

If a' = then we can set a = by redefining w, and 
from ([57)1 we find that b' / c and c'/6 are each constant, 
whence 

b c 



a 



c 
c 



where e is a constant. Writing A 
equations (|43ll46p . they become. 



(58) 

1/w in the field 



■i{A\ + e) = Kp.A^ 



kPA' 



(59) 
(60) 
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which are easily recognizable as the usual equations for 
RW spacetimes, e being the spatial curvature parameter 
usually denoted k. It is known that the most general 
solutions with zero shear, rotation and acceleration are 
the Robertson- Walker solutions (see e.g. Ellis [20l|, page 
135). Moreover from the junction condition ([26]) . i.e. 

P = 0, we have P — and consequently the fluid is a 
homogeneous collapsing dust, which has no acceleration, 
i.e. a Friedman solution. 

Hence we can state that a collapsing cylinder with 
a non-singular axis filled with shearfree irrotational 
isotropic fluid must be an RW solution and if it is 
matched to an ER solution the fluid must be dust. We 
again compare our result to the corresponding isotropic 
spherical shearfree collapse There the general solu- 
tion cannot be obtained since for the complete integration 
of the system further equations of state are required (e.g. 
an equation of state of the form P — P{^) |2l|). 

Now we show that the form for this discussed by Cocke 
is the most general one by direct coordinate transfor- 
mations. One could reach the final metric form more 
immediately by integrating (|57p and imposing regularity 
at the axis: the extra information below is that of the 
coordinate transformations. 

The RW metric can be expressed in spherical coordi- 
nates as 



ds^ = ~dt^+A^ 



dr"^ 



1 — kr'^ 



, (61) 



where A is a function only of t, fc = — 1, or 1, and the 
ranges of the coordinates are 



-oo<t<oo, 0<f<oo, 0<6'<7r. O<0<27r, 

(62) 

except when k = 1 where instead < r < 1 and r — 1 
is the antipode of the origin. To write (|6ip in cylindrical 
coordinates we make the transformation 



f = rsin0, z = /(r, 0), 



(63) 



which yields 



df = sin6'df + fcos6'd6', dz = df + d6i. (64) 



The inverse transformation is 



/ 6» df — r cos 9 dz sin d^ — / f df 

dr — . ' „ „ :r^^ w = 



svo-O ffi — r cos 9 f^f ' 



sm9 f — r cos 9 f f' 
(65) 

In order to write (|6ip transformed by (|63p we first observe 
that 



df^ 

1 — kr"^ 



+ rd9^ 



1 — kf'^ 



2f 



cos 9 
1 — kf"^ 



r sin 9 / ^ dr dz 



r^(l - fcr^ sin^ ( 



1 - kf 



dz' 



(sine f^o -rcos9 f^fY 



(66) 



Imposing the requirement that the transformation does 
not produce cross terms we must choose 



, ,_2,sin6' 
cos 9 

Substituting ([67ll back into ((66|l we obtain 

o2 a 



df2 



f^d9^ 



df^ + 



cos 



(l-fcf2)/,i 



dz' 



1 — kr'^ ■ " ^ _ j^^2 gjj^2 I 

With (IMD and ((MD we can write ^ as 
dP 



ds' ^-dt'+ A' 



1 — kf' 



+ h'dz' + f'd(j)' 



where 



h = 



cos( 



(1 - sin^ 6*) 1/2(1 - kf'y/'ff' 



(67) 



(68) 



(69) 



(70) 



Then from ^ and ^ we have 



cos ( 



/l(l-fcf2 sin2 6»)l/2(l_fcf2)l/2' 

j-^ sin 61(1 - fcf 2)1/2 



h{l - kf' sin2 61)1/2' 



Requiring that h is a function only of f we have 



nf = sm9—, nQ = r cost)—. 

dr ' dr 



(71) 
(72) 



(73) 



Differentiating ([7T|) and (|72p with respect to 9 and r re- 
spectively and using (|73p we obtain 
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f.f9 = 
f,Of = 



h[l-kf^ sin^ 0)1/2(1 _fcf2)i/2 



(1 — fcr^)sin6' rcos^ddh 



1 



h{l - kf^ sin^6')i/2 



1 - 2kf^ 



kr'^ sin 9 h df _ 

kf^ sin2 0(l-fcf2)i/2 



(l-fcf2 sin^ 6*) 1/2 



1 — kr'^ sin 9 



fl-fcr2)i/2 dh 



df 



r 



(74) 
(75) 



Equating the two expressions ([74]) and ([751) . i.e. imposing 
the integrabihty condition / ^g = f g^, and using (|63p we 
obtain 



1 d/i 
/i df 



1 — /cf 2 ' 
and integrating this we obtain 

/l = 7(l-fcf2)l/2, 



(76) 



(77) 



where 7 is an integration constant. Now substituting (j77p 
into ((69|) and rescahng z we finally have 



dP 



1 — kf' 



+ (1 - kf^)dz^ 



_ (78) 

The metric form ([75]) was obtained by Cocke @], but 
here we have proved that it is the general RW metric in 
cylindrical coordinates. 

In order to write (|78p in the coordinate system em- 
ployed in ^ we rescale t and transform f so that 



df 



(79) 



(80) 



(l-fcf2)l/2' 

and we obtain 

ds^ = A^{t){-dt^ + dr^ + g'^ dz^ + ^^2^, 
or 

ds^ = -dr^ + A'^{T){dr^ + g'^ dz^ + d^^), (81) 

with (7 — sinhr, r or sinr according as fc = — 1, or 1. 

The metric ([39]), with (gl]) and jlll), is the general RW 
metric (|80p where k — e and w = I /A, and it satisfies 
(1591) and Wl- 



VI. MATCHING FRW SPACETIME TO ER 
SPACETIME 

Matchings between a cylindrical homogeneous perfect 
fluid interior and a vacuum exterior (or vice versa) have 
been studied by Mena et al. ^] , who showed that match- 
ing to a static vacuum is impossible which is a special 
case of our result above, by Nolan and Nolan [§][3^, and 
by Tod and Mena In the last of these papers, the 
matching of ER and k — FRW metrics is given, using 
coordinates in which the exterior metric takes the form 



where R = R{f, p): the main difference from our treat- 
ment is that the coordinates (T, p) are chosen so that 
the boundary is at p = constant. Here T has replaced 
the T of the original paper to avoid confusion. Tod and 
Mena study the global structure and conclude that 
the spacetime is not asymptotically flat in the sense of 
Berger et al. [ll|, but instead has a singular Cauchy hori- 
zon. 

In , the matching of a general cylindrically symetric 
vacuum to FRW is studied: the exterior is then shown, 
as one might expect, to be of ER form. Here we shall find 
the trajectory of the boundary as an equation relating T 
and R, and then display the conditions satisfied there by 
7 and ip. Then we comment further on the results of Tod 
and Mena [13] and Nolan and Nolan ^9]. 

The well-known solutions to (l59l) and (l60l) are 



fc = 1, yl = msin^*/2, 2t = to(* - sin *), (83) 
fc =0, yl= (3V^r/2)2/3, (84) 
fc = -1, A msinh^ */2, 2r = m(sinh * - ^'J35) 

where m is a constant giving the mass density, see e.g. 
[![, equations (14.6). These are given in a form expand- 
ing as r increases from 0. Reversing the sense of r in 
([55H55|) , so that we have collapse, and introducing a con- 
stant To where the singularity occurs, we use these in 
R = A^gg' and = A{gg')' (which follow from (^5]) . 
((28)) and ([8T|) ). We can then integrate for T, and where 
necessary eliminate Vt, to get 



2a_(ro -T)= R^'^ [ - ^i?o ) {R^'^ + R^f'^ 



-i?o arcsmh 



To-T = i^oi^^/^ 
2a+(ro - T) ^ -i?i/4 ('i?i/2 + 3^ ^ (i?„ _ Ri/y 



(86) 
(87) 

/2 



-_Rn arcsin 



where Tq and Rq are integration constants, and 



s 99 



{99')" 



(89) 



ds\ ^ -e^'^'^-'f'\df^-dp^) + e^'^dz^ + e~^'''R^d(l)^, (82) evaluated on S, and thus from ([501) for e = -1 and 1 
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respectively 



s tanh ro 

a- = 



a4 



1 + tanh ro 

tanro 
1 — tan^ ro ' 



(90) 
(91) 



where r = ro, while for fc = 0, a = ro. 

On this moving boundary we know from ((24|) and 
[Ml) that we must have (noting that for all three possible 



9' 



99" 



1) 



e- = Ag', 
s 



Ar 



A^iigg'yni-i2aA,rn 

-gikjgg'y + 2Alg'^) 
Amggmi-{2aA,rr] 
9' 



(gg')Vl - (2«Ar)= 



(92) 
(93) 

(94) 

(95) 



where the functions of r have to be evaluated at rg. 

One can give a general solution of (PT|) as a sum of sep- 
arable solutions but we have not been able to determine 
the specific solution which matches to RW even in the 
simplest (fc = 0) case. 

Our condition that i? be a spatial coordinate in the 
exterior metric is violated when T is sufficiently close to 
Tq: in fact as T increases 



dR 
dT 



l4^5i?o/'(To-r) 



1/5 



(Note that this can also be expressed as 1 = |2Q;yl.T-|.) 
For larger T our coordinates, and hence our matching, 
do not apply at S. This agrees with the results of Nolan 
and Nolan who showed that such a breakdown is 
inevitable, essentially because the collapsing source al- 
ways leads to trapped cylinders, whereas the ER vacuum 
region cannot contain trapped cylinders :23|. They de- 
scribe the matching as impossible, meaning that it cannot 
be carried right up to the FRW singularity, but note that 
it can be used up to some finite time. 

One can continue the discussion of trapped cylin- 
ders by using the coordinates of Tod and Mena p^ . 
The bound on applicability of our matching corresponds 
to their conclusion that in the A: = FRW case the 
boundary becomes a marginally trapped surface when 
1 = 1 2a A, T I (in their notation, with hats added for clar- 
ity, at T = d — (4fo/3)'^) so at larger T the surface is 
trapped. Such a trapped surface is not consistent with 
"asymptotic fiatness" in the sense of Berger et al. [llj 
(see their Proposition 2.3). 

Tod and Mena further show that u — ^/Rtjj (in their 
coordinates) has a divergent derivative at the Cauchy 
horizon (essentially the past null cone of the FRW singu- 
larity) and hence conclude that this horizon is singular. 
The argument for these conclusions applies also to the 



k — ±1 cases, with changes in formulae. Tod and Mena 
infer that there is incoming gravitational radiation: how- 
ever, this does not seem to be supported by the calcula- 
tions in the following section. 



VII. ENERGY, SUPERENERGY, RADIATION 
AND BOUNDARY CONDITIONS 



Part of our motivation was a search for an exact cylin- 
drical solution for interior and exterior enabling one to 
study exactly how gravitational radiation arises. Our 
ansatze for the interior turned out to allow only FRW, 
which we would expect to be non-radiative in any defini- 
tion. Thus we expect the exterior to also be non-radiative 
(or possibly to show waves coming in from infinity and 
totally refiected at the boundary with the interior). 

A cylindrically symmetric spacetime cannot be asymp- 
totically fiat, due to the behaviour in directions parallel 
to the axis, so one cannot employ the usual global defi- 
nition of radiation for isolated bodies due to Bondi et al. 
p^ . Moreover, it appears from Proposition 2.3 of Berger 
et al. [ll| that simple cylindrical solutions with collaps- 
ing cores could not even be "asymptotically flat" in their 
modified sense, because one would expect trapped cylin- 
drical surfaces to arise in any collapse which does not 
halt or reverse, and such surfaces are not compatible with 
"asymptotic fiatness" . 

So for more detailed study we need some definition of 
radiation, or energy, or energy density, other than the one 
from asymptotic fiatness. This must be (quasi-)local, at 
least in z, to avoid the problem that integration for z 
from — oo to oo would obviously give an infinite answer 
for any non-zero energy density. 

There is also the possibility of energy being transferred 
to or from any given region by transport along the axis. 
Bondi [l| showed that there is no conserved mass per unit 
length as a result of 'intangible' gravitational induction 
arising from axial motion (and distinct from the work 
done by axial pressure). 

It is well known that in relativity there cannot be any 
covariant local definition of energy, as this would violate 
the equivalence principle. There are a number of quasi- 
local definitions, using either integrals over surfaces, or 
pseudo-tensors integrated over volumes (which in gen- 
eral also reduce to surface integrals, provided the inte- 
rior volumes do not have discontinuities or singularities of 
the pseudo-tensor). Clavering j24i] has shown that none 
of the latter agree with the various quasi-local surface 
integrals discussed by Szabados 0, and none of them 
are satisfactory, for a variety of reasons, the best being 
MoUer's definition. We therefore do not calculate the 
pseudo-tensorial energies. 

In his study of definitions of standing waves, Stephani 
(26| considered cylindrical systems. His results suggest 
that for the ER solutions Thome's "C-energy" ^] is the 
least unsatisfactory. (A recent indication of the unsatis- 
factoriness of C-energy, even in vacuum, has been given 
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in [27| . where it is shown that it can be non- vanishing 
in Minkowski space.) We consider, following Chiba |28| 
and Hayward 29], the modified C-energy defined in the 
"Note added in proof" on pages B256-257 of Thome's pa- 
per. Taking the generic cylindrically symmetric metric in 
the form 

ds^ = -e^^-'-^\dT^-dr^) + e^^dz^ + e~^^R^d^^, (96) 
where R — R{r, t), this energy is 

i(l-(i?'2_7?2)e-27) (97) 

where the prime and dot refer to differentiation with re- 
spect to r and t respectively. 

Hayward has shown, in an elegant formulation, that 
one can define an invariant tensor (in his notation, 
Qab) such that the sum of this and the usual energy- 
momentum tensor Tab is conserved: may then be inter- 
preted as the energy-momentum of gravitational waves. 
One can derive the conservation in a simple way from the 
fact that for the metric ([M]) . when matter is present, the 
equations generalize to 

i[(i?'2 - i?2)e-27]' ^ -RR'{nToo + i''^ + (98) 

+ RR{kToi+2'iP'iP), 
\ [(i?'2 - i?2)e-27]- ^ -^RR\nToi + (99) 

+ RRinTn + i^'^ + V-^) 

(cf. equations (A41), (A42), of [29]): the result is then 
just the integrability condition for the left sides, the rel- 
evant terms of Q being the terms in ^ on the right sides. 
That these terms are invariantly defined (provided there 
are no more translational Killing vectors) arises because 
the two Killing vectors, 9^ and 9^, and their lengths, are 
uniquely defined up to rescaling of z, so is fixed up to 
a constant and i? up to a constant factor. 

Hayward's discussion also brings out the fact that for 
the form (|96p . if fi and p are the divergences of the in- 
coming and outgoing null normals to timelike cylinders 
of symmetry, then 

P^^=^e'^^-^\R'^~R^) (100) 

with an obvious relation to (|97|) . 

This is also related to one of the best known quasi- 
local energy definitions by an integral on a surface, the 
Hawking mass. For a closed surface S with surface area 
element dS this mass is 

To make a closed cylindrical surface we would have to 
add (e.g.) surfaces z — constant at some finite z values, 
but jip has the same values for all z on such surfaces and 
so we could ignore them for very large cylindrical surfaces 



and think of E as giving a Hawking mass per unit length 
in z. 

Unfortunately this does not lead to a satisfactory ac- 
count of energy lost or gained by the collapsing dust. The 
Darmois conditions imply that E is continuous at E (one 
can check this by direct calculation, but it is easy to un- 
derstand because there is no surface layer and hence no 
immediate change in geodesic deviations at S, so pp is 
continuous, and the continuity of and gzz then shows 
the same is true for E.) Calculating on the dust side of 
E, using ([SO)) , we have 

(As one expects, [1^, this is ©(r^) as r ^ 0.) At a fixed 
r, only the / = A\ term changes, and it increases 
as the dust collapses. This is consistent with Cocke's 
calculation Q of (unmodified) C-energy flux, normal to 
cylinders of constant i?, in our notation, in the Einstein- 
Rosen region, where he found an inward flux. 

It would, however, be physically very odd if this had to 
be interpreted as the exterior giving energy to the dust, 
since the behaviour of the dust is exactly the same as 
in a uniform universe, where a cylinder cannot easily be 
thought of as taking energy from the rest of the universe. 
Probably the result is better considered as another indi- 
cation of the unsatisfactoriness of C-energy. 

Other local characterizations of the presence of radia- 
tion are given by the decomposition of the Bel-Robinson 
tensor. The behaviour of this tensor in an ER spacetime 
was considered in [s^], where the obvious unit vector in 
the form (PO)) . i.e. the one parallel to d/dT, was used to 
define the decomposition. In that reference it was shown 
that for a pulse of radiation, well behind the front of the 
pulse, there is an incoming flux of superenergy, which 
seems to agree with Cocke's result. However, to define 
states of intrinsic radiation one has to consider all possi- 
ble timelike vectors, and we now show, using the timelike 
vector parallel to the boundary, that our solutions do not 
have intrinsic radiation there according to this definition. 

Bel [111, HI] defined the tensor 

j.abcd ^ R<^ecfj^bd ^ *^aecf*^J,d 

+R*aecfj^*^b^d ^ *j^*aecf*j^*^b^d^ ^^^^^ 

which in the vacuum case is referred to as the Bel- 
Robinson tensor. Here the star operation is the usual 
Hodge dual (see e.g. @, chapter 3). Bonilla and Sen- 
ovilla (33I pointed out that with the usual decomposition 
in terms of the Weyl and Ricci tensors, Cabcd and Rab, 
i.e. 

Rabcd — Cabcd + Eabcd + Gabcdj 

(102) 
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where 

Eabcd = \{gacSbd + ObdSac ~ 9adSbc ~ ffbc^ad) , (103) 
Gabcd = YiRiSacgbd — gaddbc) = -^R-Qabcdi (104) 
Sab = Rab — jRgab, R = R"'ai (105) 

the Bel tensor can be written as 

rpabcd ^aecf ^ b ^ , ^ ^aecf ^ ^ b ^ 

^R,^acbd^(jadbc^^^j^abcd^ (--^Qg^ 

6 

where the matter contribution is 

Mabcd = S^'S''" 

+ lS'"'S''eg'='' + IS'^'S'^eg"^ - 2S<''g''^'^''Sf 

+i5'=/5e/(g''^5'' + g^'g"' - g^'g'") (i07) 

+ ^(25 + 2g"V^ - 5"V')- 

One can decompose the Bel tensor in a 3+1 formahsm 
[s^ l relative to a unit timelike vector n°, giving, among 
the parts, the superenergy W, the superPoynting vector 
and the tensor Qabc defined by 

W = T"''""'nanbn,nd, P'' = T^''''nanbn,h\, 
Qbcd ^ ^T''''f3nah\h''fh'^g, (108) 

where hab = gab + narn,. Bel [l^l defined a state of in- 
trinsic gravitational radiation (at a point p) to be one in 
which P° 7^ for any choice of (at p). Garcia-Parrado 
Gomez-Lobo js^l similarly defines an intrinsic superen- 
ergy radiative state to be one where Q'^''^ ^ for any 
choice of In vacuum, 

Qcdb = h,dPb ~ 2{BdaEcf - B.aEdsW^ . (109) 

where, as usual, Eab = Cacbdn^n'^, Bab = *Cacbdn^n'^, 
and -qabc = rjabcdn'^, Vabcd being the usual volume 4- form 
density. 

We now consider how these quantities behave at the 
(timelike) boundary between regions of spacetime. It 
is simplest to describe this using an orthonormal tetrad 
{Ca, a = 0...3} chosen such that the timelike unit vector 
Co lies on the boundary surface (and will be used as 
in p08p ) and ei is the normal to the surface. From the 
equations (80) of Mars and Senovilla 35] we can straight- 
forwardly show that if the Darmois junction conditions 
are satisfied then the following combinations of Riemann 



tensor components are continuous: 

Ei2, Ri2, Ei3, i?i3, Rqi, (110) 

Bii, B22 (and hence B33), P23, (HI) 

B13 + 2^02, B12 — 2R03, E23 — 2R23 (112) 

—£'11 ^ ^R + \{R22 + R33), E22 + \R + \{Roo — R22), 

E33 + ^R + ^{Roo — ^33)- (113) 
Note that one consequence is (f26| . the continuity of Tn. 

From (jllip and (|112p we see that (assuming we align 62 
and 63 with Sa and K^) then with the energy- momentum 
form assumed in ([1]), Bab is continuous at the boundary. 
Moreover, the reflection symmetries in and z imply that 
P2 = -P3 = (and P12 = P13 = 0) in both interior and 
exterior (note that since reversing one axis also reverses 
the orientation and hence the sign of the dual, it does 
not follow that P23 — 0). 

From this, it is easy to see that for an FRW interior, 
which is conformally flat, Bab = on both sides of the 
boundary, so from (|109p Pa = and Qabc = there in the 
frame defined above (i.e. with n" = Cq and under Dar- 
mois boundary conditions). Hence at the boundary the 
exterior does not have intrinsic gravitational radiation 
or intrinsic superenergy radiation. Thus it is reasonable 
to conclude that any possible radiation in the exterior 
spacetime is not produced by the source. This agrees 
with our expectation that such an interior should not 
radiate or absorb radiation. Note that we have not ex- 
cluded the possibility of total reflection at the boundary 
if at a reflecting surface these indicators (P° and Qcdb) 
would show no intrinsic radiation, and it is therefore pos- 
sible that there could be non-zero incoming and/or out- 
going radiation in the exterior. The conclusion that no 
radiation crosses the boundary seems in conflict with the 
discussion given by Tod and Mena [l3|. 
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